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Reduction of thermal noise in dielectric mirror coatings is a key issue for the sensitivity improve- 
ment in second and third generation interferometric gravitational wave detectors. Replacing an end 
mirror of the interferometer by an anti-resonant cavity (a so-called Khalili cavity) has been proposed 
to realize the reduction of the overall thermal noise level. In this article we show that the use of a 
Khalili etalon, which requires less hardware than a Khalili cavity, yields still a significant reduction 
of thermal noise. We identify the optimum distribution of coating layers on the front and rear sur- 
faces of the etalon and compare the total noise budget with a conventional mirror. In addition we 
briefly discuss advantages and disadvantages of the Khalili etalon compared with the Khalili cavity 
in terms of technical aspects, such as interferometric length control and thermal lensing. 
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I. INTRODUCTION 

The sensitivities of second-generation (Advanced 
LIGO, Advanced VIRGO, GEO-HF, and LCGT) and 
third-generation (Einstein Telescope) interferometric 
gravitational wave detectors will be partly limited by 
thermal fluctuations in the mirrors [1-5]. 

The pioneering articles on this issue were dedicated 
to the investigation of the mirror substrate fluctua- 
tion: Brownian thermal noise [6-8] and thermo-clastic 
noise [9]. Fundamental thermal motion (Brownian mo- 
tion) of material atoms or molecules causes Brownian 
noise. Fundamental thermodynamic fluctuations of tem- 
perature lead to thermo-elastic noise through the ma- 
terial's thermal expansion. Similarly thermo-refractive 
noise [10, 11] is caused by temperature fluctuations lead- 
ing to fluctuations of the refractive index and therefore 
fluctuations of the optical path length inside the mate- 
rial. These results were obtained for the model of an 
infinite test mass, i.e. the mirror was considered to be an 
elastic layer with infinite width and finite thickness. All 
of these results were generalized for a finite-size mirror 
model [8, 12]. 

Very soon the importance of mirror coating thermal 
noise was realized as its parameters may differ consid- 
erably from the mirror substrate parameters. Despite 
its low thickness, the very high loss angle of the mirror 
coating materials (usually Si02 and Ta205) makes the 
coating Brownian noise the most significant one among 
all kinds of mirror thermal noise [13-15]. The thermo- 
elastic noise of the coatings only has a small contribu- 
tion in the total noise budget [16, 17]. Later, Kimble [18] 
proposed the idea of thermal noise compensation which 
was explored carefully in [19, 20] for a particular case of 



thermo-elastic and thermo-refractive noise. 

The coating Brownian noise is still one of the main 
contributions to the noise spectra of gravitational wave 
observatories [13-15]. One of the most promising ap- 
proaches aimed to decrease its level was offered by Khalili 
[21] who proposed to replace the end mirror in the inter- 
ferometer arm with a short Fabry-Perot cavity tuned to 
anti-resonance (see center panel of Figure 1). In prac- 
tice most light is reflected from only a few first layers 
(farest from substrate) and all others (located closer to 
substrate) only rcfiect a small part of light. However, 
since the thermal fluctuations are proportional to the to- 
tal thickness of the coating and the inner layers of coating 
are the main contribution to the phase fluctuations of the 
reflected light, the transmittance of each mirror can be 
higher to realize the same reflectivity of the system as a 
compound end mirror. The total thickness of coatings in 
such Khalili cavity mirrors is the same as the thickness 
of the conventional mirror, while the Brownian noise of 
the end mirror of a Khalili cavity is significantly reduced 
[21] because the thickness fluctuations of the second end 
mirror coating (EEM in Figure 1) do not influence the 
fluctuations of the input mirror coating (lEM in Figure 
1). Moreover, using a rigidly controlled Khalili cavity 
allows a reduction of coating Brownian noise [22]. 

One of the main problems in the Khalili cavity is to 
establish a low-noise control of the mirror positions (see 
detailed explanation in Sec. III). A potentially easier way 
is to use a Khalili etalon (KE) instead of a Khalili cavity 
(KG) or a simple conventional mirror (GVI). The idea is 
to use a single mirror but to split the coating into two 
parts (see right hand panel of Figure (EEM in Figure 1): 
the front coating (on the front substrate surface) features 
just a few layers and the rear coating (on the rear sub- 
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Figure 1: Simplified schematic of an Advanced LIGO interferometer with conventional end mirrors, featuring all the coating 
layers on their front (left), with Khalili cavities as end mirrors (center) and with Khalili etalons, featuring only a few coating 
layers on the front surface and the majority of the coating layers on the rear surface (right). 



strate surface) consists of the rest of the required coating 
layers. 

The purpose of this article is to develop an idea of 
the Khalili etalon [23] , to calculate the total mirror ther- 
mal noise arising in a KE and in a CM, and to compare 
them. We investigate the idea of using a KE in the Ein- 
stein Telescope (ET) and Advanced LIGO (aLIGO). In 
Sec. II we describe the mirror parameter optimization 
procedure, namely the optimal number of layer pairs in 
the front coating. In Sec. II D we describe the details 
of the thermal noises arising in the KE and CM calcula- 
tions. Section III is dedicated to the problem of thermal 
lensing which is much more important in a KE than in 
a CM. In Sec. IV we discuss the obtained results and 
draw the conclusions. Finally, some calculation details 
are provided in the appendices A-B. 



II. COATING OPTIMIZATION 

The main idea of using a KE is to reduce the mirror's 
total thermal noise without reducing its reflectivity. By 
total thermal noise spectral density we mean the sum of 
the Brownian, thermo-elastic and thermo-refractive noise 
spectral densities. Coating Brownian noise is caused 
mostly by the fluctuations of the entire coating thick- 
ness. It would then seem evident that Brownian noise 
be lower when the front coating contains less layers and 
hence the lowest noise be achieved for the coating totally 
displaced to the rear mirror surface. This is in princi- 
ple true but at the same time some other noises, such as 
substrate thermo-refractive noise, rise dramatically caus- 
ing the total noise level to rise also. Moreover, the less 
layers one puts onto the front coating, the higher will 
be the absorption in the substrate. So there has to be 
an optimum of how to best distribute the coating layers 
between the front and back surfaces in order to obtain 
minimal total thermal noise and not too much of absorp- 



tion in substrate. The aim of this section is to find this 
optimum configuration. 



A. Thermal Noise calculation technique 

The only way to find the optimal number of front coat- 
ing layers, A'^i, is to compare the thermal noise for ev- 
ery A^i. This requires the calculation of the different 
noise contributions as functions of the front coating lay- 
ers number Ni. The most basic principles we used are: 
(i) the total number of Ta205 and Si02 layers, N, is 
fixed, i.e. we used the coating structure planned for both 
ET and aLIGO and modified it to fit the double coating 
paradigm: 20 Ta205 and 18 Si02 quarter-wave layers 
plus the substrate (it is considered as an ordinary but 
"slightly" thicker coating layer) and plus two caps con- 
sisting of a half-wavelength Si02 layer (for the CM it 
would have been 20 Ta205 layers and 19 Si02 layers plus 
one cap); (ii) a quarter- wavelength Ta205 layer and a 
quarter-wavelength Si02 layer are alternately coated on 
the front or rear surface so that there are always an odd 
number of front coating layers (TVi = 1, 3, 5 etc.) and also 
an odd number of layers of the rear coating N2 = N — Ni 
(37, 35, 33 etc.). Please note that the substrate and 
caps arc not included in these numbers; (iii) the num- 
ber of layers of the front coating A^i is the argument 
and the total thermal noise driven mirror displacement 

•^totai i^ t^^*^ function of it; (iv) we consider only Brow- 
nian, thermo-elastic and thermo-refractive noises (being 
the most significant contributions), and (v) wc used the 
mirror of a finite-size cylinder, the model of which has 
been developed in Ref. [8, 12, 15, 23], and calculated 
all noises numerically using the fluctuation-dissipation 
theorem (FDT) [7, 24, 25] as it is briefly described in 
Sees. IIAl-IIC. The optimal number of front coating 
layers appeared to be A^i = 3, i.e. 2 Ta205 layers and 1 
Si02 layer plus a cap in the front coating and 18 Ta205 



layers and 17 Si02 layers plus a cap in the rear coat- 
ing. With the technical feasibility taken into account 
(see Sec. II C), however, it turns out that the system 
with A^i = 5 layers on the front surface (i.e. 3 Ta205 
layers and 2 Si02 layers plus a cap on the front mirror 
surface and 17 Ta205 layers and 16 Si02 layers plus a 
cap on the rear surface) will be better and we analyze 
the system with TVi = 5 in detail. In this case the mirror 
thermal noise does not reach its minimum but it is only 
about 5 % higher. 



1. Brownian noise 

The total coating thermal noise of the etalon will be 
the sum of noise on the front surface and noise on the 
back surface: 



Sx = eiSxi + 626x2 



(1) 



Here 6x1 is the displacement of the front surface of the 
mirror and 6x2 the displacement of the boundary sur- 
face between the rear surface of the mirror substrate and 
the coating on it. Considering the KE as a Fabry-Perot 
cavity consisting of two mirrors with amplitude reflectiv- 
ities Ri (front coating) and R2 (rear coating) tuned to 
anti-resonance, one can calculate the coefficients ei and 
£2 (see details in Rcf. [23]): 



ei 



£2 = 



i?i [1 + (1 + n,)RiR2 + i?i] + i?2(l - Us) 



(l+i?li?2)2 



n,R2{l-Rl) 
{I + R1R2Y 



(2) 



Here Ug is the substrate refractive index. Note that Ri 
and R2 are functions of the number of front and rear 
coating layers. In particular, we have the following for- 
mulas for i?i and R2 as functions of the number of the 
front coating layers A^i (A''2 = A^ — TVi is the number of 
the rear coating layers): 
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where rii and 712 arc Ta2 05 and Si02 coating layers 
refractive indices. 

Hence, in order to calculate spectral density S{uj) of 
the displacement Sx caused by thermal noise using the 
FDT, one has to apply the forces 

eiFoe*"* and £2^06*"* 

to the front and rear coatings correspondingly and to cal- 
culate the total dissipated power [7, 24, 25]. For the cal- 
culation of the spectral density S'b(w) of Brownian coat- 
ing noise the dissipated power may be calculated through 



the elastic energy f/g ' stored in each fc-th layer (of the 
front or rear coating): 

+E^J'M';'^+E^^Mf)dr, (4) 

where Eij and Ty are the strain and stress tensor com- 
ponents (only the non-zero components are shown in the 
formula above) , R is the mirror radius and h^ is the thick- 
ness of the fc-th layer. The components Eij and Tij are 
calculated as it is described in detail in [23]. Then the 
Brownian noise spectral density may be evaluated as fol- 
lows: 
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where ks is Boltzmann's constant, T is the absolute tem- 
perature and (j>k is the loss angle describing structural 
losses in the k-th layer. The sum is taken over all layers, 
i.e. TV = 38 is the number of layers without the substrate 
and the caps. The total number of summands is there- 
fore A^ -I- 3 = 41; A^ = 38 layers in the front and rear 
coatings, plus 2 layer-caps and 1 layer-substrate. So the 
index fc = refers to the front coating cap, the indices 
A: = 1 -^ A^i refer to the front coating quarter wavelength 
(QWL) layers, the index A: = iVi + 1 refers to the sub- 
strate, the indices k^Ni + 2^N + l refer to the rear 
coating QWL layers and the index k = N + 2 refers to 
the rear coating cap. There are A^ -I- 2 = 40 summands; 
the ones with k ^ Ni are to be considered for coating 
Brownian noise and the one with index fc = A^i -I- 1 is to 
be considered for substrate Brownian noise: 
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Here (f>Ni+i = 0s is the loss angle of the substrate, while 
01 and 02 represent the loss angels of the Ta205 and 
Si02 layers, respectively. The values are presented in 
Table II. 



2. Thermo-elastic noise 

The thermo-elastic (TE) noise calculations for the sub- 
strate and for the coating is similar. In order to calculate 
the dissipated power one should calculate the diagonal 
components Ejj of the strain tensor for each layer (in- 
cluding the substrate and the caps) and take the trace: 



eW=i^^)+i?^+i?W. 



Then one may find the power dissipated through the TE 
mechanism [12, 25]: 
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where k^ is the thermal conductivity, Ck is the thermal 
capacity per unit volume, Yk is the Young's modulus, 
Vk is the Poisson's ration, a^ is the thermal expansion 
coefficient, pk is the density, and the index k denotes 
number of the layer. Therefore, the TE noise spectral 
density will simply be: 
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fe=0 



Similar to the Brownian noise calculations, the sum- 
mands with the indices fc ^ TVi + 1 are relevant for the 
coating TE noise while the one with the index fc == A^i + 1 
needs to be considered for the TE noise of the substrate: 
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3. Thermo-refractive noise 

TR noise originates from thermodynamic fluctuations 
of the temperature ST in the substrate, producing phase 
fluctuations of the reflected wave phase via the tem- 
perature dependence of the substrate's refraction index 
Us- Likewise, the phase fluctuations may be recalculated 
into effective fluctuations of mirror surface displacement 
6x = —e2Psh5T where the coefficient £2 introduced in 
(2), characterizes the light amplitude circulating inside 
the substrate and /3s = dus/dT is the thermo-optic coef- 
fleient of the substrate. 

We calculate the thermo-refractive (TR) noise in the 
substrate using the model of an infinitely large plane in 
the transverse directions with thickness of h. The spec- 
tral density of the temperature fluctuations in this model 
is shown in [26] see Eq. (E8): 
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where w is the radius of the light spot (intensity decreases 
with distance r from center as ^ e~'^^ /™ ) and the pa- 
rameters with subscript s refer to the substrate. The TR 
noise spectral density for the substrate (recalculated to 
displacement) becomes 
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Benthem and Levin have pointed out that some correc- 
tions should be applied to this formula. This corrections 
are based on the account of the fact that light inside the 
arm froms the standing wave and not a traveling wave. 
We can rewrite Eq. (2) of Ref. [11] in a simpler form with 
only the normal incidence and the circular beam being 
considered: 
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In addition we have to consider the TR noise present 
in the coatings [10]. For its estimate we use the following 
formula 
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(12) 
(13) 



where A is the wavelength of light in vacuum, P^s is the 
averaged thermo-optic coefficient of the entire coating, 
and /?! = dni/dT and /?2 = dn2/dT are the thermo- 
optic coefficients of Ta205 and Si02 layers, respectively. 
This formula is based on the assumption that only the 
flrst few layers contribute considerably to the thermo- 
refractive loss mechanism. It is obtained for a mirror 
with an infinite radial dimension and a finite height. This 
model is valid with good accuracy for CM. However, for 
KE we use the same formula as an order-of-magnitude 
estimation. 



B. Optimization results 

In this subsection we present the results of our opti- 
mization process. Using the proposed parameters we ob- 
tained numerical estimates of all noise sources discussed 
above for ET and aLIGO. All geometrical design param- 
eters for these interferometers are presented in Table I. 
The physical constants and material parameters are sum- 
marized in Table II. 

First of all, we analyze the spectral density of the dis- 
placement noise 5'jo^jji for a Khalili etalon (KE) as a 
function of the number of front layers A'^i. This noise 
analysis considers the sum of the noise sources listed in 
Sec. II A: Brownian, TE and TR noises which are divided 
into a coating and a substrate contribution each. The KE 
total thermal noise is then compared to the results for a 



conventional mirror (CM) S'lojaj 
G which is defined as: 
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using a gain parameter 



(14) 
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This gain G has to be maximized in order to enhance 
the detector sensitivity. In Fig. 2 we plot the gain G 



as a function of the number of front coating layers A'^i 
for both, ET and aLIGO. Please recall that the number 
of rear coating layers iV2 = A^ — iVi is constrained by 
the total number of coatings iV = 38. One can see that 
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Figure 2: Gain G (see (14)) in ET (blue stars) and aLIGO (red 
stars) as a function of the number of front coating layers A''i . 
The number of rear coating layers N2 results from the total 
number of layers A'" = 38 minus the number of front layers A'^i . 
Parameters used for this calculation are presented in Tables 
I and II. For ET a maximum gain of G = 1.72 appears for 
Ni = 3 front coating layers. For aLIGO the maximum gain 
of G = 2.16 is also obtained for A''i = 3. 

for both detectors, ET and aLIGO, the gain is obviously 
maximized for the case of A'^i =3 front coating layers. 
Then the maximum gain appears to be G = 1.72 for ET 
and G = 2.16 for aLIGO. 



C. Absorption 

Another important parameter to be taken into account 
is the substrate absorption *abs- It describes the portion 
of light energy which is absorbed in the substrate with 
respect to the light incident to the mirror. The value of 
^abs is accessible via the light power circulating inside 
the substrate. We know that this light power inside the 
etalon will be a factor of 

1-Rl 



{l + RiR^f 



lower than the light power incident to the mirror. Finally, 
the value vE'abs may be evaluated using the absorption 
coefficient of the substrate material rj (in ppm per cm) 
and the substrate thickness h (in cm): 
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(I + R1R2) 



(15) 



Note the factor of 2 in front of the substrate thickness, 
which occurs as the light passes the substrate twice: once 



forward and once backward. Also keep in mind that the 
reflectivities i?i and R2 are functions of the number of 
front and rear coating layers — see formulas (3). 




V> „15 20 .25 30 3! 

JN umber oi tront coating layers A'l 



Figure 3: Substrate absorption *abs in ET (blue stars) and 
aLIGO (red stars) as a function of the number of front coating 
layers TVi . The number of rear coating layers N2 results from 
the total number of layers A'' = 38 minus the number of front 
layers A'^i. Parameters used for this calculation are presented 
in Tables I and II. 

As Fig. 3 illustrates ^abs decreases exponentially with 
an increasing number of front coating layers. For the 
optimum number of front coating layers A^i = 3 the sub- 
strate absorptions in ET and aLIGO equal 2.67 ppm and 
1.78 ppm, respectively. It seems reasonable to assume 
that a loss coefficient of 1 ppm is admissible. In this case 
we have to choose A^i = 5. Indeed, using formula (3) 
with A'^i = 5 and coating parameters listed in Tables I 
and II we obtain: 
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(16) 



It means that for aLIGO (circulating power Wq 



0.8 MW) the absorbed power is about W, 
as for ET (Wo = 3 MW) — W^ 



(ET) 
abs 



(AL) 

abs 



0.72 W, 



4.1 W. 



Consequently, everywhere below in this article we as- 
sume the number of front coating layers to be A'^i =5. 
This choice allows the gain to be G = 1.67 for ET and 
G = 2.03 for aLIGO. 



D. Thermal noise of the Khalili etalon 

In Fig. 4 and 5 we present the thermal noise spec- 
trum of a KE including different noise sources for ET 
and aLIGO, respectively. For numerical estimates we use 
the parameter data listed in Tables I and II. One clearly 
realizes that Brownian thermal noise dominates the mir- 
ror thermal noise at almost all the frequency range from 
1 Hz to 10 kHz. 



Table I: Parameters used for the numerical calculation, w is 
the radius of the laser beam (intensity decreases with distance 
r from the center as ~ e"^"" '^ ), R and h are radius and 
thickness of the cylindrical mirror, Wo is the power circulating 
in the arm cavities of the interferometer. 



Parameter 


Einstein Telescope 


Advanced LIGO 


w, m 


0.12 


0.06[17] 


R, m 


0.31 


0.17 


h, m 


0.30 


0.20 


Wo, MW 


3 


0.8[27] 


^(CM) 


39+cap 


39+cap 


iV(KE) 


38 + 2 caps 


38 + 2 caps 



Table II: Parameters used for the numerical calculations. T 
is the temperature, A is the optical wavelength in vacuum, 
rj denotes the optical losses per unit length in the substrate 
and n is the refractive index. The remaining parameters are 
explained in Sec. IIA1-IIA3. 



Parameter 


substrate 


TaaOs layer 


Si02 layer 


T, K 




300 


[17] 


A, m 


1 


.064 X lO"'' 


[17] 


V, ppm/m 


25 


- 


- 


n[17] 


1.45 


2.035 


1.45 


/?, 1/K [17, 19] 


8 X 10-" 


1.4 X 10"^ 


8 X 10-" 


Q, 1/K[17] 


5.1 X 10"^ 


3.6 X 10~" 


5.1 X lO-'^ 


p, kg/m3[17] 


2202 


6850 


2202 


Y, Pa 


72 X 10^ [17] 


140 X 10^ [28] 


72 X 10'''[17] 


u 


0.17[17] 


0.23[28] 


0.17[17] 


K, W/K m[17] 


1.38 


33 


1.38 


C, J/K kg[19] 


746 


306 


746 


<^ 


4x 10-1" [29] 


2 X 10-*[30] 


4x 10-*[31] 



Sepcrately, we also present the numerical results for all 
noise sources at a single frequency of / = 100 Hz in Ta- 
ble III. We choose this frequency as a round number lo- 
cated in the frequency range with the highest sensitivity. 
This frequency value has already been used previously in 
this article in Sec. II B for numerical estimates. 

Brownian noise is the main object of our investigations 
as it dominates the sensitivies of both detectors (ET and 
aLIGO) almost in the whole detection band. It can be 
calculated accurately using the model developed in [23]. 
Shortly, a calculation method is presented in Sec. II A 1. 
A further inspection of the spectral noise sensitivity plots 
reveals substrate Brownian noise to be the second impor- 
tant noise process. Thus, Brownian noise dominates over 
TE and TR noise absolutely and is the main factor lim- 
iting both interferometers sensitivities in the frequency 
domain near 100 Hz. 

Taking the same noise sources into account as for our 
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Figure 4: ET noise spectral densities for KE with a front 
coating of A'^i =5 layers (plus a cap) and a rear coating 
of N2 = 33 layers (plus a cap). Parameters used for this 
calculation are presented in Table I and Table II. 



KE investigation (see Sec. IIA1-IIA3) we have calcu- 
lated the numerical values for the spectral noise density 
of a corresponding CM. For clarity we do not present 
the single contributions but only the total mirror ther- 
mal noise in this section. Using thermal noise values for 



,(KE) 



(CM). 



KE and CM [S^^^ ' and S^^^ ') and the definition of the 
gain parameter G (14) we compare KE and CM to esti- 
mate the benefit of using a KE. For the spectral densities 
calculated in Table HI we arrive at a gain of 

Get = 1.67, for ET parameters, (17) 

Gal = 2.03, for aLIGO parameters. (18) 

The sligtly larger gain value for aLIGO parameters may 
be qualitatively explained by the mirror geometry. So 
the ET mirror is more sensitive to membrane deforma- 
tions than the aLIGO mirror. This property allows to 
'transfer' Brownian fluctuations from the rear coating to 
the front surface more effectively. Indeed, the geometri- 
cal factor g = 2R/h (the fraction of diameter to thickness 
of mirror) for ET is larger than for aLIGO: 
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1.7, 
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Semi-qualitative consideration 



In this subsection we would like to present a way to 
simply estimate the gain. For an order of magnitude 
estimate we may approximate the total thermal noise 
by Brownian coating noise that prevails at all frequency 
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Figure 5: aLIGO noises spectral densities for KE with a front 
coating of A^i =5 layers (plus a cap) and a rear coating 
of A''2 = 33 layers (plus a cap). Parameters used for this 
calculation are presented in Table I and Table II. 



Table III: Numerical values for all noise sources (Brownian, 
TE and TR noise of mirror substrate and coating) and total 
mirror thermal noise at 100 Hz for ET and aLIGO. Parame- 
ters used for this calculation are presented in Tables I and II. 
The front coating of the KE consists of A''i = 5 layers (plus a 
cap) while the rear coating consists of A'^2 = 33 layers (plus a 
cap). 



noise spectral density 


Einstein Telescope 


Advanced LIGO 


KhaliUEtalon (KE): 
coating Brown., -^ 
substrate Brown., -§= 

' vHz 

coating TR, ^ 
substrate TR, ^= 

VHz 

substrate TE, ^= 

VHz 

coating TE, ^ 


1.70 X 10"^^ 
6.06 X 10-^2 
3.58 X 10-^2 
1.60 X 10-22 
7.98 X 10-23 
1.79 X 10"2* 


2.63 X 10-21 
8.55 X 10-22 
7.16 X 10-22 
3.38 X 10-22 
2.14 X 10-22 
6.63 X 10-2-1 


KE total, ^= 

VHz 


1.85 X 10-21 


2.89 X 10-21 


Conv. Mirror (CM): 
coating Brown., -^= 
substrate Brown., -S= 

VHz 

coating TR, ^ 
substrate TE, ^ 
coating TE, ^^ 


2.99 X 10-21 
6.53 X 10-22 
3.58 X 10-22 

8.51 X 10-23 

4.52 X 10-2'' 


5.75 X 10-21 
9.32 X 10-22 
7.16 X 10-22 
2.31 X 10-22 
1.78 X 10-23 


CM total, ^= 

' vHz 


3.09 X 10-21 


5.88 X 10-21 



ranges as we have seen. Moreover, the thickness of Si02 
layers and Ta205 layers differs only about 40 % while the 
Ta205 loss angle is 5 times higher than the loss angle of 
Si02 . Therefore, we may very roughly approximate the 
total thermal noise with the sum of Ta2 05 coating layers 
Brownian noise (recall that all of them are uncorrelated) . 
In a first approximation one could assume that the front 
coating is responsible for the main contribution to the 
noise level. Thus, the total thermal noise level should be 
proportional to the number of front coating Ta205 layers 
only. For a CM this number is Qcm = 20 and for a KE 
~ Qke"' = ^- ^^^ ratio of the Q values should represent 
the gain of a KE: 



Gn 



Q front 
KE 



f-- 



(20) 



The estimated gain is larger compared to the accurate 
results (17). It may be explained by the fact that we do 
not account for elastic coupling (through substrate) be- 
tween rear coating layers motion and front coating layers 
motion, i.e. the displacement of the front coating due to 
a deformation of the rear coating layers. One could say, 
the rear coating layers motion is 'transferred' to the front 
coating through the substrate. This coupling is moder- 
ated by the clastic properties of the latter. 

We introduce a transfer ratio p to account for this elas- 
tic coupling. The variable p ranges from for a Khalili 
cavity (KC) to 1 for a CM or 'zero'-thickness substrate. 
We can calculate p using the simple model of a cylindrical 
mirror whose front and rear surface are covered by equal 
layers (same thickness and same elastic parameters). Let 
us apply a single force at the front surface and keep the 
rear surface free of forces. One can calculate the elas- 
tic energies in the front layer C/front and in the rear layer 
t^rear- Obviously, the transfer ratio p may be calculated 
as the ratio of both energies. This estimate gives: 



C/r< 



P 



Ufr, 



PET ^ 0.23, PAL ^ 0.086 (21) 



The ratio pal for aLIGO is smaller than for ET. Again 
this behaviour can be explained by the different geometry 
factors g (see estimates (19)). 

Now instead of Eq. (20) we can state a more accurate 
formula for the gain estimate taking into account the 
elastic coupling p of the Q^e ^ 1''' ^'^^^ coating layers: 



Ga 



g,(ET) 
app 



/CM 



/^front 
"^KE 



-pxQ™-' 



1.71, 



^(AL) 
app 



2.12. 



(22a) 
(22b) 



We see that the approximated gain values coincide with 
the accurate values (17) within an accuracy of about 5 %. 
Real gains in ET and aLIGO are lower than the expected 
approximated values (22) because of other noise sources 
that were omitted here (Brownian substrate and Brown- 
ian coating of the Si02 layers, substrate and coating TE 
and TR noise). 



Note that the elastic couphng docs not take place in 
a KC where both coatings arc mechanically separated 
by vacuum. Consequently for both detectors, ET and 
aLIGO, the usage of a KC instead of a CM is expected 
to show a gain value of G « 2.6 (as (20)). 



E. Potential sensitivity improvements for future 
GW interferometers 



In this section we quantitatively analyse the overall 
sensitivity improvement potentially achievable by replac- 
ing the conventional end mirrors by KE in aLIGO and 
ET. 
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Figure 6: Potential overall sensitivity improvement of aLIGO 
for the use of KE, compared to conventional end mirrors: 
The binary neutron star inspiral range increases by 15 % from 
196.2 to 225.0 Mpc. This corresponds to a potential increase 
in the detected event rate for binary neutron star inspirals of 
slightly above 50 %. 

Figure 6 shows the potential sensitivity improvement 
of aLIGO for the use of KC as end mirrors. The sensitiv- 
ity curves have been created using the GWINC software 
[32] and for a signal recycling configuration that is opti- 
mised for the detection of binary neutron star inspirals 
(see configuration 2 in [33]). Only the two main noise 
contributions arc shown: quantum noise (black trace) 
and coating Brownian noise (sum of all test masses) (red 
trace), as well as the total noise (blue traces). Please 
note that all other relevant noise sources have been in- 
cluded in the calculations of the total noise traces, but 
have been omitted from the plot for clarity. The dashed 
lines indicate the strain levels for the standard aLIGO 
design, while the solid lines show the potentially reduced 
noise levels originating from the application of KE as end 
test masses, as described in this article. The main dif- 
ference between these two scenarios originates from the 
reduction of coating Brownian noise by a factor 2.18, as 
described by the values in the right hand column of Table 
III. Please note that thermal noise contributions from the 
input mirrors stay identical for the two scenarios. The 



Figure 7: Potential overall sensitivity improvement of the ET 
high frequency detector (as described in [34]) for the use of 
KEs, compared to conventional end mirrors; The binary neu- 
tron star inspiral range increases by 15 %, corresponding to a 
potential increase in the detected event rate for binary neu- 
tron star inspirals of 50 %. 



corresponding increase in the binary inspiral range (1.4 
solar masses, SNR of 8, averaged sky location) is about 
15 % and therefore yields an relative increase of the bi- 
nary neutron star inspiral event rate of about 50 %. 

Figure 7 shows the sensitivity improvement of a poten- 
tial ET high frequency detector as described in [34] for 
the replacement of the conventional end mirrors by KEs. 
Following the values given in Table III we considered a 
flat coating Brownian noise reduction factor of 1.76 for 
the end test masses, while again we assumed the ther- 
mal noise of the input test masses to stay constant. This 
yields an overall reduction of the total thermal noise of 
all test masses of about 25 % and an increase in the ob- 
servatory sensitivity of up to 20 % in the most sensitive 
frequency band between 50 and 400 Hz. We find an in- 
crease in the binary neutron star inspiral range of 15 % 
from 1593 to 1833 Mpc. This corresponds to an increase 
in the binary neutron star inspiral event rate of about 
50%. 



III. TECHNICAL FEASABILITY OF KHALILI 
ETALONS FOR FUTURE GW OBSERVATORIES 

A. Required Hardware 

Figure 1 shows the simplified schematics of an aLIGO 
or ET interferometer with different end mirror configu- 
rations. Replacing the conventional end mirrors by KCs 
would have a significant impact on the required hard- 
ware. Instead of a single end mirror suspended from a 
single seismic isolation system per end mirror, in the case 
of the KC two mirrors with two full seismic isolation sys- 
tems are required at the end of each arm cavity. 

This means that with KCs there are six (2x IM, 2x 



lEM, 2x EEM) instead of four optical elements (2x IM, 
2 EM) , which require the maximal seismic isolation. The 
concept of the KE allows us to still significantly reduce 
the thermal noise contribution of the end mirrors, while 
being compatible with the already available seismic isola- 
tion systems. Therefore, it is in principle possible to up- 
grade a 2nd or 3rd generation gravitational wave detector 
by replacing conventional end mirrors by KEs without 
altering or extending the vacuum systems and scimsic 
isolation systems. 



B. Interferometric Sensing and Control 

In addition to the reduced hardware requirements, the 
main advantage of the KEs with respect to KCs is the 
potential simplification of several aspects related to the 
interferometric sensing and control. Upgrading aLIGO or 
an ET interferometer from its standard configuration to 
employ KCs increases the length degrees of freedom of the 
main interferometer from five (DARM (differential arm 
length), MICH (Michelson cavity length), SRCL (signal 
recycling cavity length), CARM (common arm length), 
PRCL (power recycling cavity arm length)) to a total of 
seven. 

It is worth mentioning that the additional two degrees 
of freedom actually have a very strong coupling to the dif- 
ferential arm length channel of the interferometer. For 
the example of the coating distribution discussed in this 
article, the length of the KC needs to be stabilized with 
an accuracy of only a factor 10 less than what is required 
for the main arm cavities. That means the length of the 
KC needs to be orders of magnitude more stable than for 
example the differential arm length of the central Michel- 
son interferometer. 

In order to achieve this demanding stability of the KC 
one has to make use of highly dedicated readout and 
control schemes. Special care needs to be taken to avoid 
introducing potential control noise at the low frequency 
end, which could potentially spoil the overall sensitivity 
of the gravitational wave detector. 

Substituting the KC, consisting of two individual mir- 
rors potentially encountering independent driven motion 
(for example seismic), by the proposed KE would ensure 
that both relevant mirror surfaces would be rigidly cou- 
pled via the etalon substrate. Therefore, the length of 
the KE would be much less susceptible to seismic dis- 
turbances or gravity gradient noise, as compared to the 
length of KC. Also in terms of potential control noise the 
KE is advantageous over the KC. In case of the KC the 
mirror positions would have to be controlled by means of 
coil magnet actuators or electro-static actuators, which 
can potentially introduce feedback noise at frequencies 
within the detection band of the gravitational wave de- 
tector. In contrast the length of the KE can be locked 
by controlling the etalon's substrate temperature (using 
the temperature dependency of the index of refraction). 
Since the etalon substrate acts as a thermal low pass, the 



etalon length will be extremely constant for all frequen- 
cies within the detection band of the interferometer. 

However, not only the length sensing and control is 
highly demanding in case of a KC, but also the align- 
ment sensing and control. Again the key point here is 
to find a high signal to noise error-signal and then ap- 
plying low noise feedback systems to keep the mirrors of 
the KC aligned in pitch and yaw. As the KC would be 
rather short compared to the main arm cavities, the KCs 
would unfortunately feature a high mode degeneracy, i.e. 
it would not only be resonant for the desired TEMqo 
mode, but also for higher order modes, which would fur- 
ther increase the alignment requirements. Using a KE 
would potentially allow us to transfer the alignment con- 
trol problem from the detector operation to the manu- 
facturing process of the etalon. If it would be possible 
to manufacture an KE with sufficiently parallel front and 
back surface, we would not need to actively control the 
relative alignment of the KE surfaces during operation. 
The two parameters that would be most relevant are the 
relative curvature mismatch of the etalon front and rear 
surfaces as well as the parallelism of the two surfaces. 
As we have shown in [35] the curvature mismatch is the 
dominating factor for the etalon's performance. 



C. Thermal Lensing 

In this section we will compare the thermal lensing [39] 
of the KC configuration to the one of a KE. In the case 
of the KC we have the following absorption processes: 
(i) lEM front coating, (ii) lEM substrate, (iii) lEM anti- 
reficx coating on its rear surface and (iv) front coating 
of EEM. In the case of the proposed KE the situation 
is pretty similar apart from process (iii), which docs not 
exist. 

In the following we will show by means of FEM that 
the actual thermal lensing induced into the KE is of the 
same order, but slightly smaller than in the case of the 
KC. 

The FEM used here treats the mirror as a substrate. 
The coatings and the laser beam are included as heat 
sources. For the reflective coatings we assume an absorp- 
tion of 0.5 ppm, and 1 ppm for the anti-reflective coating 
in the KC. The FEM assumes an emissivity e = 0.93, an 
ambient temperature of 300 K and uses the parameters 
from Table I and II for the values of aLIGO. After com- 
puting the temperature and displacements of the flnite 
elements, the optical path difference (OPD) is derived. 
For the OPD wc included the temperature dependence 
of the refractive index (which is the dominant thermal 
lensing effect in fused silica) and the expansion, while we 
omitted the elasto-optic effect. We also did not include 
surface to surface radiation in the KC. This would make 
the thermal lens worse, and is therefore safe to exclude 
in order to make a conservative comparison between KC 
and KE. Fig. 8 shows the temperature distribution in 
a KC and a KE for the aLIGO parameters with Ni=5, 
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Figure 8: The temperature distribution in the mirrors of a KC (left) and a KE (right) for aLIGO parameters (with Ni=5) as 
calculated via FEM. Please note the different color scales for the lEM and the EEM in the KC. 



OPD cross section for a Khalili cavity 



OPD cross sections in a Khalili etalon 




-10 -5 5 10 

Distance from ttie optical axis in cm 



-0.5 



' 






/ \ 


dn/dT 


/ \. 


dV/dT 


/ v 


total 


/ \ 


thin lens fit 


/ \ 


f=2004m 


/ 


K 


A K- 



-15 -10 -5 5 10 15 

Distance from tfie optical axis in cm 



Figure 9: The optical path delay (OPD) is computed in the FEM simulation, in order to determine the focal length of the 
thermal lens in the KC (left) and the KE (right) for aLIGO parameters. Shown are the total OPD, the OPD caused by the 
thermo-refractive effect (named dn/dT) and the OPD due to the expansion of the substrate {dV/dT). The plots show also fits 
of the OPD for an ideal thin lens. 



while Fig. 9 presents the corresponding OPD for a sin- 
gle pass due to the thernio-optic effect and expansion of 
the substrate, as well as a fit of an OPD that would be 
caused by an ideal thin lens. The fits are least square 
fits, weighted by the beam intensity. For the aLIGO pa- 
rameters with Ni=5, the thermal Icnsing in the KC can 
be described by a thermal lens with a focal length of 
/ = 1944 m, while the OPD in the KE can be fitted by 
a thermal lens with a focal length of / = 2004 m. The 
respective values for ET arc a focal length of 1797 m for 
the KC and 1838 m for KEs. It follows that the induced 
thermal lensing in the KC and KE is of similar strength. 

As we have shown above, the thermal lensing for the 



KE is slightly weaker than for the KC. 

As one can see from the magnitude of the induced 
thermal lensing, the compensation of this effect will be 
extremely challenging in both cases. Potential ways of 
mitigating the thermal lensing could include innovative 
approaches such as radiative cooling [38] or pre-shaped 
mirror or etalon substrates, which feature the wrong cur- 
vature, when being cold, but develop the correct 'shape' 
when operated at the designed optical power. The com- 
pensation of thermal lensing has turned out to be more 
challenging than anticipated in the first generation grav- 
itational wave detectors. Only, the practical experience 
that will be collected with the advanced detectors will 
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allow us to realistically judge the feasibility of KCs as 
well as KEs. However, the main purpose of the ther- 
mal lensing analysis presented here was to show that the 
thermal lens will not be worse, but slightly better for the 
proposed KE compared to a KC. 



IV. CONCLUSION 

In this article we have investigated the main thermal 
noise sources arising in the mirrors of the two next gener- 
ation gravitational wave detectors: Advanced LIGO and 
Einstein Telescope. The thermal noise sources include 
Brownian, thermo-elastic and thermo-refractive noise of 
the mirror coatings and the mirror substrate, among 
which the Brownian coating noise is the largest. We ap- 
plied our model developed in [23] to study the idea of 
the Khalili etalon to decrease the coating thermal noise 
and to improve the sensitivity. The optimum KE con- 
figuration minimizing the total therixial noise level was 
found to be with 2 Ta205 layers and 1 Si02 layer plus 
a cap in the front coating and with 18 Ta205 layers and 
17 Si02 layers plus a cap in the rear coating. However, 
since the substrate absorption in ET with such a configu- 
ration is 8.00 W, and that in aLIGO is 1.42 W, our choice 
is not to use the optimal but a slightly different coating 
distribution with A^i = 5, i.e. 3 Ta205 and 2 Si02 lay- 
ers plus a cap on the front surface and 17 Ta205 and 
16 Si02 layers plus a cap on the rear surface. The ab- 
sorbed power in the substrate with such a configuration is 
4.1 W for ET and 0.72 W for aLIGO. Such an absorption 
is around 1 ppm which seems to be reasonable price for 
the thermal noise enhancement. The total noise spectral 
density of ET and aLIGO can be improved by the factors 
of 1.67 and 2.03, respectively, compared with the cases 
of conventional end mirrors . Moreover, we have checked 
our numerical calculations with a very simple qualitative 
consideration designed to make an order of magnitude 
estimation. This estimation shows an agreement of bet- 
ter than 5 percent with the exact numerical calculations. 
A use of KEs instead of conventional end mirrors would 
improve the detection rate of the binary neutron star in- 
spirals with the future gravitational wave observatories 
by about 50%. 

We also discussed the feasibility of the Khalili etalon 
compared with that of the Khalili cavity. The KE is more 
advantageous in terms of the hardware requirements. We 
also compared the thermal lensing effects in the KE and 
in the KC and found that the former is slightly better 
for not having the anti-reflective coatings that KC con- 
tains on the rear surface of the front mirror. In fact, the 
thermal lensing problem in either case is quite severe and 
we should explore a way to compensate the lensing effect 
without imposing excess noise. 

In this paper we assumed that the light is reflected 
on the outer surface of each coating without taking into 
account reflections from inner layers. A more accurate 
analysis shown in Ref. [36] gives a value of coating Brow- 



nian noise slightly lower (about 10%). 

Thermoelastic and thermo-refractive noises originate 
from a thermodynamical fluctuation of the temperature 
and the correlation of the two noises can be non-trivial 
with a certain set of parameters [19, 20]. In this pa- 
per, the correlation was ignored and we treated the two 
noises individually, which is not a problem as one of them 
is much lower than the other in the case of KE (see Ta- 
ble HI). It should be noted, however, that the optimal 
KE configuration could be determined in such a way that 
thermoelastic noise and thermo-refractive noise be neg- 
atively compensated if the mechanical loss angles of the 
coating materials were 10 times lower than the current 
values. 
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Appendix A: Coefficients ei and £2 calculation 

Let us consider KE as a Fabry-Pcrot interferometer 
with two mirrors (namely two reflective coatings) with 
amplitude transmi ttances Ti and T2, and am plitude re- 
flectivities i?i = ^1 - T2 and i?2 = \/l - T|. The mir- 
rors are separated by a medium with a refractive index 
Us, and the mean distance between the mirrors is L. Op- 
tical losses are equal to zero. The fluctuations of the 
coordinates of the front and rear mirrors are reprenseted 
by X and y, respectively. The probe beam A is incident on 
the front mirror (coating) and is partially reflected. We 
are interested now in the reflected beam B. The weight 
coefficients ei and 62 represent how much the fluctuations 
X and y contribute to the reflected beam B, respectively. 

For a short cavity we can use a quasi-static approxi- 
mation ~ it means that the motion of the mirrors are suf- 
ficiently slow compared with the relaxation rate of the 
cavity. We assume that the optical path between the 
mirrors is fixed to a quarter wavelength, i.e. e^^"^'^ = i. 

We can consider the cavity as a generalized mirror. Ob- 
viously, the reflectivity of the generalized mirror depends 
on the fluctuations x and y. However, for the reflected 
beam we have to include the motion of the generalized 
mirror, that is, the common-mode motion of x and y. 
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The reflected beam shall be described as: 

R2^l - Ri 



B 



2ikx 



A 



1 - RiR2'&l ' 



^r, 



2ikns {y—x) 



(Al) 



We have already taken into account the fact that e**^"''^ = 
i (cavity is tuned in the anti-resonance) and t?o describes 
the variation of cavity length due to the mirror fluctua- 
tion. 

The fluctuations x and y arc small enough compared 
with the cavity length L so that we may expand (Al) 
into series over x and y. Keeping the linear terms only, 
we get 
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(A3) 

(A4) 



{I + R1R2) 
which have been introduced in Eq. (2). 



Appendix B: Coating reflectivities Ri and R2 
calculations 



Let us first consider a CM with N altering layers of 
Ta205 and Si02 with the refractive indices ni and n2, 
respectively, and the substrate with the refractive index 

Us- 

A multilayer coating consisting of layers with refractive 
indices nt and lengths k is described by the same formu- 
las for the transmission line consisting of ports with wave 
resistances l/rii and the distances li [37]. It is convenient 
to describe the transmission line with impedances Zi and 
reflectivities Ri . Impedance Zi of the i-th layer can sub- 
stitute the total impedance of all the layers between this 
layer and the substrate, which does not affect the other 
layers. It is convenient to start the calculation from the 
boundary of the substrate and the A'^-th layer, and Zq 
will be the equivalent impedance of all the mirror and 
Rq will be the reflectivity of the entire system. 

There is a recurrent formula for impedances and re- 
flectivities of neighboring layers (neighboring ports of the 
transmission line): 



where 0, 
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Akriili 
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Ri — 
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(Bl) 

(B2) 



is the phase shift in the i-th layer. 



Using (Bl) and (B2) one may easily get the recursive 
formula 
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(B3) 



The substrate is considered as an infinite half-space, 
so its impedance is given by Zjv = 1/njv+i = l/f^s 
and hence its reflectivity is given by Rn = (1 — 
nAr/n7v+i)/(l + nAr/njv+i) = {ns-nN)/ins+nN). Then 
Eqs. (B3) and (B2) yields each Zi and Ri. We are in- 
terested only in Rq. The thickness of each layer in the 
high-reflective coating is a quarter- wavelength (QWL), 



I.e. 
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Then (B3) becomes: 
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Using (B4) one may easily get the chain (remember 
that N is odd and refractive indices alter so that rii = 
713 = ■■■ = nN = ni, n2 = 714 = • • • = "w-i = "2): 
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(B5b) 
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and thus the total coating reflectivity is: 
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(B6) 



The cap does not change the impedance Zq. The length 
of the cap Ic is a half- wavelength (HWL) so that e''^^"^" = 
— 1. Using (B3) one may see that the impedance of the 
system with the cap is exactly the same as that without 
the cap: Zc-i = Zc. At last, Eq. (B6) for Rq is valid 
with or without the cap. 

Now we are interested in refiectivities of the Khalili 
etalon coatings (that should be used as refiectivities Ri 
and i?2 in a Fabry-Perot interferometer used in Ap- 
pendix A). For both coatings we may use the same for- 
mula (B6) but with different number of layers A'^i and 
N2 = N — Nl and different 'border' refractive indices Ug 
and uq: 

For the front coating both the substrate and the vac- 
uum play their original roles i.e. the rear and front infi- 
nite half-spaces respectively. Thus, Hs = n2, no = 1 and 
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N^Ni: 
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Ha I "2 
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For the rear coating the vacuum plays the role of 
the substrate (the rear infinite half-space) and the sub- 
strate plays the role of the vacuum (the front infinite half 



space). Thus, rig — I, Uq = n2 and N — N2: 

N2-1 



i?2 



"2 ( "2 



1 



12 ( 22 ) 



No-l' 



(B8) 



If we use the fact that N2 = N ~ Ni these formulas 
(B7-B8) will coincide with formulas (3). 
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